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We present a search for a narrow-width heavy resonance decaying into top quark pairs (X ti) 
in pp collisions at =1.96 TeV using approximately 0.9 fb"'^ of data collected with the DO detector 
at the Fermilab Tevatron Collider. This analysis considers it candidate events in the lepton plus 
jets channel with at least one identified b jet and uses the ti invariant mass distribution to search 
for evidence of resonant production. We find no evidence for a narrow resonance X decaying to tt. 
Therefore, we set upper limits on ax ■ B{X ^ tt) for different hypothesized resonance masses using 
a Bayesian approach. For a Topcolor-assisted technicolor model, the existence of a leptophobic Z' 
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boson with mass Mz' < 700 GeV and width Tz' = 0.012Mz> can be excluded at the 95% C.L. 



PACS numbers: 14.65.Ha, 14.70.Pw 

INTRODUCTION 

The top quark has by far the largest mass of all the 
known fermions. Unknown heavy resonances may play 
a role in the production of top quark pairs (tt) and add 
a resonant part to the standard model (SM) production 
mechanism mediated by the strong interaction. Such res- 
onant production is possible for massive Z-like bosons in 
extended gauge theories [1], Kaluza-Klein states of the 
giuon or Z boson [2, 3], axigluons [4], Topcolor [5], and 
other theories beyond the SM. Independent of the exact 
model, resonant production of top quark pairs could be 
visible in the reconstructed ti invariant mass distribution. 

In this Letter, we present a search for a narrow-width 
heavy resonance X decaying into ti. We consider the 
lepton+jets (^+jets, where £ = e ot ^) final state. The 
event signature is one isolated electron or muon with high 
momentum transverse to the beam axis (pr), large trans- 
verse energy imbalance {^t) due to the undetected neu- 
trino, and at least four jets, two of which result from the 
hadronization of b quarks. The analyzed dataset corre- 
sponds to an integrated luminosity of 913±56 pb"^ in the 
e-|-jets channel and 871 ± 53 pb""*^ in the /i+jets channel, 
collected with the DO detector between August 2002 and 
December 2005. The analysis uses events with at least 
three reconstructed jets. Backgrounds from light-quarks 
are further reduced by identifying b jets. After b tag- 
ging, the dominant physics background for a resonance 
signal is non-resonant SM it production. Smaller contri- 
butions arise from the direct production of W bosons in 
association with jets (VF-|-jets), as well as instrumental 
background originating from multijet processes with jets 
faking isolated leptons. The search for resonant produc- 
tion in the ti invariant mass distribution is performed 
using Bayesian statistics to compare SM and resonant 
production to the observed mass distribution. 

Previous searches performed by the CDF and DO col- 
laborations in Run I found no evidence for a fl reso- 
nance [6, 7]. In these studies, a Topcolor model was 
used as a reference to quote mass limits. According to 
this model [5], a large; top quark mass can be generated 
through the formation of a dynamical ti condensate, Z' , 
due to a new strong gauge force with large coupling to 
the third generation of fermions. In one particular model, 
Topcolor-assisted technicolor [8], the Z' boson has large 
couplings only to the first and third generation of quarks 
and has no significant couplings to leptons. Limits ob- 
tained on ax ■ B{X ^ ti) are used to set a lower bound 
on the mass of such a leptophobic Z' boson. In Run I 
CDF found Mz' > 480 GeV with 106 pb~^ of data [6], 
and DO obtained Mz' > 560 GeV using 130 pb"^ [7], 



both at the 95% C.L. and for a resonance with width 
Tz' = O.OUMz'. 

DO DETECTOR 

The DO detector [9] has a central-tracking system 
consisting of a silicon microstrip tracker and a central 
fiber tracker, both located within a 2 T superconducting 
solenoidal magnet, with designs optimized for tracking 
and vertexing at pseudorapidities |?7| < 3 and \r]\ < 2.5, 
respectively. The pseudorapidity, 77, is defined with re- 
spect to the beam axis. Central and forward preshower 
detectors are positioned just outside of the superconduct- 
ing coil. A liquid-argon and uranium calorimeter has a 
central section (CC) covering pseudorapidities \r]\ < 1.1, 
and two end calorimeters (EC) that extend coverage to 
r/l « 4.2, with all three housed in separate cryostats [10]. 
An outer muon system covering |ry| < 2 consists of a layer 
of tracking detectors and scintillation trigger counters in 
front of 1.8 T iron toroids, followed by two similar layers 
after the toroids [11]. Luminosity is measured using plas- 
tic scintillator arrays placed in front of the EC cryostats. 
The three-level trigger and data acquisition systems are 
designed to accommodate the high luminosities of Run II 
and record events of interest at up to about 100 Hz. 

EVENT SELECTION 

To select top quark pair candidates in the e-|-jets and 
^-l-jets decay channels, triggers that required a jet and 
an electron or muon are used. The event selection re- 
quires either an isolated electron with pt > 20 GeV 
and \ri\ < 1.1, or an isolated muon with px > 20 GeV 
and |?7| < 2.0. No additional isolated leptons with 
Pt > 15 GeV are allowed in the event. Details of the 
lepton identification and isolation criteria are described 
in [12, 13]. We require to exceed 20 GeV (25 GeV) 
for the e-|-jets (^+jets) channel. Jets are defined using 
a cone algorithm [14] with radius T^cone = 0.5, where 
^conc = y (A(^)2 + (Ay)2, (f) is the azimuthal angle, and 
y the rapidity. The selected events must contain three or 
more jets with pt > 20 GeV and |y| < 2.5. At least one 
of the jets is required to have pT > 40 GeV. Events with 
mismeasured lepton momentum are rejected by requiring 
the to be acollinear with the lepton direction in the 
transverse plane: A0(e,^T) > 2.2 - 0.045 GeV~^ and 
A(j){i2,^T) > 2.1 - 0.035 GeV"^.^T [15]. 

To improve the signal-to-background ratio, at least one 
jet is required to be identified as a 6 jet. The tagging al- 
gorithm uses the impact parameters of tracks matched 
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to a given jet and information on vertex mass, the de- 
cay length significance, and the number of participating 
tracks for any reconstructed secondary vertex within the 
cone of the given jet. The information is combined in 
a neural network to obtain the output variable, NN^, 
which tends towards one for b jets and towards zero for 
light quark jets [16]. In this analysis we consider jets to 
be 5-tagged if NN^ > 0.65 which corresponds to a tag- 
ging efficiency for b jets of about 55% with a tagging rate 
for light quark jets of less than 1%. 

We independently analyze events with three and four 
or more jets and separate singly tagged and doubly 
tagged events, since the channels have different signal- 
to-background ratios and systematic uncertainties. 

SIGNAL AND BACKGROUND MODELING 

Simulated events are used to determine selection ef- 
ficiencies for the resonant ti production signal and for 
background sources except those in which instrumental 
effects give fake leptons and in multijet production 
events. Samples of resonant tt production are generated 
with PYTHIA [17] for ten different choices of the resonance 
mass Mx between 350 GeV and 1 TcV. In all cases, the 
width of the resonance is set to Fx = 0.012M;!s:- This 
qualifies the X boson as a narrow resonance since its 
width is smaller than the estimated mass resolution of 
the DO detector of 5-10%. The generated resonance is 
forced to decay into tt. 

Standard model tt and diboson backgroimds (WW, 
WZ, and ZZ) are generated with pythia [17]. Single 
top quark production is generated using the COMPHEP 
generator [18]. A top quark mass of 175 GeV is used for 
both resonant and SM top production processes. VF-l-jets 
and Z+jets events are generated using alpgen [19] to 
model the hard interaction and pythia for parton show- 
ering, hadronization and hadron decays. To avoid double 
counting between the hard matrix element and the par- 
ton shower, the MLM jet-matching algorithm is used [20]. 
The CTEQ6L1 parton distribution functions [21, 22] are 
used for all samples. The generated events are processed 
through the full GEANT3-bascd [23] simulation of the DO 
detector and the same reconstruction program as used 
for data. 

The SM tt, single top quark, diboson, and Z-|-jets 
backgrounds are estimated completely from Monte Carlo 
(MC) simulation, to obtain the total acceptance as well 
as the shape of the reconstructed tl invariant mass dis- 
tribution. Trigger inefhciencies and differences between 
data and MC lepton and jet identification efficiencies are 
accounted for by weighting the simulated events [15]. 
Jet 6-tagging probabilities are measured in data and 
parametrized as functions of pr and 77. They are used 
to weight each simulated event according to its event 
fo-tagging probability. Finally, the expected yields are 



TABLE I: Event yields for the expected SM background and 
for data. The uncertainties are statistical and systematic. 





3 jets 


> 4 jets 


tt 


167.4 


160.5 


VF-l-jets 


118.2 


24.1 


Other MC 


34.8 


9.8 


Multijet 


31.3 


7.4 


Total background 


351.7±29.3 


201.8 ±29.0 


Data 


370 


237 



normalized to the SM theoretical prediction. A tt pro- 
duction of atj = 6.77 ± 0.60 pb for mt = 175 GeV [24] 
is used. Z-|-jets, single top quark and diboson samples 
are normalized to their next-to-leading-order cross sec- 
tions [25 27]. 

The VF-|-jets background is estimated from a combina- 
tion of data and MC information. The expected number 
of VF+jets events in the 6-tagged sample is computed 
as the product of the estimated number of I^-l-jets be- 
fore b tagging and the expected event 6-tagging proba- 
bility. The former is obtained from the observed number 
of events with real leptons in data, computed using the 
matrix method [12], and then subtracting the expected 
contribution from other SM production processes. The b- 
tagging probability is obtained by combining the W-|-jets 
flavor fractions estimated from MC with the event b- 
tagging probability, estimated from b tag rate functions. 
The shape of the reconstructed invariant mass distribu- 
tion is obtained from the MC simulation. 

The multijet background is completely determined 
from data. The total number of expected events is esti- 
mated by applying the matrix method to the each of the 
6-tagged subsamples. The shape is derived from events 
with leptons failing the isolation requirements. A sum- 
mary of the prediction for the different background con- 
tributions in the combined £+jcts channels, along with 
the observed number of events in data, is given in Ta- 
ble I. Systematic uncertainties are discussed below. 

RECONSTRUCTION OF THE tt INVARIANT 
MASS DISTRIBUTION 

The ti invariant mass is reconstructed from the four- 
momenta of up to the four highest pT jets, the lepton mo- 
mentum, and the neutrino momentum. The latter is ob- 
tained from the transverse missing energy and a W-mass 
constraint. The neutrino transverse momentum is iden- 
tified with the missing transverse momentum, given by 
and its direction. The neutrino momentum along the 
beam direction, p'^, is estimated by solving the equation 
— (P^ where p^ {p") is the lepton (neutrino) 

four momentum. If there are two solutions, the one with 
the smaller is taken; if no solution exists, p'^ is set to 
zero. This method gives better sensitivity for high mass 
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FIG. 1: Shape comparison of the expected tt invariant mass distributions for SM top quark pair production (histogram) and 
resonant production from narrow-width resonances of mass Mx = 450, 650, and 1000 GeV, for (a) 3 jets events and (b) > 4 jets 
events. 



resonances than a previously applied constrained kine- 
matic fit technique [7], while only slightly reducing the 
sensitivity for lower resonance masses. Moreover, this di- 
rect reconstruction allows the inclusion of data with fewer 
than four jets in the case that some jets are merged, fur- 
ther increasing the sensitivity. The expected ti invariant 
mass distributions for three different resonance masses 
are compared to the SM expectation in Fig. 1. 



SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties can be classified as those 
affecting only normalization and those affecting the shape 
of any of the signal or background invariant mass distri- 
butions. The systematic uncertainties affecting only the 
normalization include the theoretical uncertainty on the 
SM prediction for a^j (9%), the uncertainty on the in- 
tegrated luminosity (6.1%) [28], and the uncertainty on 
the lepton identification efficiencies. 

The systematic uncertainties affecting the shape of the 
invariant mass distribution as well as the normalization 
are studied in signal and background samples. These in- 
clude uncertainties on the jet energy calibration, jet re- 
construction efficiency, and 6-tagging parameterizations 
for b, c and light jets. The effect due to the top quark 
mass uncertainty is computed by changing TOj in the sim- 
ulation of tt to 165 GeV and 185 GeV, normalized to their 
corresponding theoretical cross sections. The effect is 
scaled to correspond to a top quark mass uncertainty of 
±5 GeV. The difference in the ti acceptance due to the 
top quark mass variation is also included in the system- 
atic uncertainty. 

The fraction of heavy flavor in the ly-t-jets background 
is measured in control samples, and a corresponding un- 
certainty on the VF-|-jets flavor composition is used. Also 



the uncertainties on the 6-fragmentation and the uncer- 
tainties of the efficiencies used in the matrix method are 
taken into account. 

Table II gives a summary of the relative systematic un- 
certainties on the total SM background normalization for 
the combined £-|-jets channels. The effect of the different 
systematic uncertainties on the shape of the tt invariant 
mass distribution cannot be inferred from this table, but 
is included in the analysis. 

RESULT 

After all selection cuts, 319 events remain in the e-|-jets 
channel and 288 events in the fi+jets channel. The sums 
of all SM and multijet instrumental backgrounds are 



TABLE II: The relative systematic uncertainties on the over- 
all normalization of the SM background and for a resonance 
mass of Mx = 650 GeV, with at least one b-tagged jet. The 
uncertainties shown are symmetrized. The actual asymmet- 
ric uncertainties and the effect of shape-changing systematic 
errors are used in the limit setting. 





SM processes 


Resonance 


Source 


(backgrounds) 


Mx = 650 GeV 




3 jets 


> 4 jets 


3 jets 


> 4 jets 


Jet energy calibration 


±1.0% 


±5.8% 


±3.7% 


±5.5% 


Jet energy resolution 


±0.2% 


< 0.1% 


±1.2% 


±0.2% 


Jet identification 


±0.6% 


±2.0% 


±0.6% 


±1.6% 


attimt = 175 GeV) 


±3.1% 


±5.9% 






Top quark mass 


±5.2% 


±6.9% 






b tagging 


±3.1% 


±4.9% 


±3.9% 


±3.6% 


b fragmentation 


±0.3% 


±0.4% 


±0.6% 


±0.6% 


Vy-|-jets (heavy flavor) 


±2.5% 


±0.9% 






Multijet lepton fake rate 


±0.3% 


< 0.1% 






Selection efflciencies 


±3.1% 


±5.3% 


±3.6% 


±3.6% 


Luminosity 


±2.6% 


±4.2% 


±6.1% 


±6.1% 
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FIG. 2: Expected and observed ti invariant mass distribution for the combined (a) 1 + 3 jets and (b) ^" + 4 or more jets channels, 
with at least one identified b jet. Errors shown on the data points are statistical. Superimposed as white area is the expected 
signal for a Topcolor-assisted technicolor Z' boson with Mz' ~ 650 GeV. 



303 ±22 and 251 ±19 events, respectively. Invariant mass 
distributions are computed for events with exactly one b 
tag and for events with more than one b tag. Addition- 
ally, the distributions are separated into 3 jets and > 4 
jets samples. The measured invariant mass distributions 
and corresponding background estimations are shown in 
Fig. 2 for the 3 jets and > 4 jets samples. 

Finding no significant deviation from the SM expecta- 
tion, we apply a Bayesian approach to calculate 95% C.L. 
upper limits on ax'B{X ^ti) for hypothesized values of 
Mx between 350 and 1000 GeV. A Poisson distribution 
is assumed for the number of observed events in each 
bin, and flat prior probabilities are taken for the signal 
cross section times branching fraction. The prior for the 
combined signal acceptance and background yields is a 
multivariate Gaussian with uncertainties and correlations 
described by a covariance matrix [29]. 

The expected and observed 95% C.L. upper hmits on 
ax-B{X ^ti) as a function of Mx, after combining the 1 
and 2 6-tag samples and the 3 and > 4 jets samples, are 
summarized in Table III and displayed in Fig. 3. This 
figure also includes the predicted ax-B(X ^ti) for a lep- 
tophobic Z' boson with Tz' = 0.012Mz' computed using 
CTEQ6L1 parton distribution functions. The compari- 
son of the observed cross section limits with the Z' boson 
prediction excludes Mz' < 700 GeV at the 95% C.L. Due 
to a small excess of data over expectation (of no more 
than 1.5f7 significance) for invariant masses in the range 
between 600 and 700 GeV, the observed limits do not 
reach the expected limit for a Z' boson of 780 GeV. 



TABLE III: Expected and observed limits for ax-B{X -*tt) 
at the 95% C.L. when combining all channels and taking all 
systematic uncertainties into account. 

'Mx [GeV] Exp, limit [pb] Obs. limit ^pEf 



350 


2.08 


3.19 


400 


2.09 


2.32 


450 


1.59 


1.59 


500 


1.24 


0.99 


550 


0.94 


0.80 


600 


0.68 


0.79 


650 


0.55 


0.87 


750 


0.36 


0.66 


850 


0.28 


0.49 


1000 


0.22 


0.36 



5f 




§00 400 500 600 700 800 900 1000 1100 



Mx[GeV] 

FIG. 3: Expected and observed 95% C.L. upper limits on 
ax-B{X —>ti) compared with the predicted Topcolor-assisted 
technicolor cross section for a Z' boson with a width of F^/ = 
0.012Mz' as a function of resonance mass Mx- The shaded 
band gives the ±1 sigma uncertainty in the SM expected limit. 
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CONCLUSION 

A search for a narrow-width heavy resonance decay- 
ing to tt in the ^+jcts final states has been performed 
using data corresponding to an integrated luminosity of 
about 0.9 fb~^, collected with the DO detector at the 
Tevatron collide^'. By analyzing the reconstruc;ted tt in- 
variant mass distribution and using a Bayesian method, 
model independent upper limits on ax ■B{X ^ti) have 
been obtained for different hypothesized masses of a 
narrow-width heavy resonance decaying into ti. Within 
a Topcolor-assisted technicolor model, the existence of a 
leptophobic Z' boson with Mz' < 700 GeV and width 
Tz' = 0m2Mz> is excluded at the 95% C.L. 
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